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Abstract: 

We present details of a calculation of the cross section for hadronic top-antitop produc- 
tion in next-to-leading order (NLO) QCD, including the decays of the top and antitop into 
bottom quarks and leptons. This calculation is based on matrix elements for 1/^^'^ ij,~i'fjhh 
production and includes all non-resonant diagrams, interferences, and off-shell effects of 
the top quarks. Such contributions are formally suppressed by the top-quark width and 
turn out to be small in the inclusive cross section. However, they can be strongly enhanced 
in exclusive observables that play an important role in Higgs and new-physics searches. 
Also non-resonant and off-shell effects due to the finite W-boson width are investigated 
in detail, but their impact is much smaller than naively expected. We also introduce a 
matching approach to improve NLO calculations involving intermediate unstable parti- 
cles. Using a fixed QCD scale leads to perturbative instabilities in the high-energy tails of 
distributions, but an appropriate dynamical scale stabilises NLO predictions. Numerical 
results for the total cross section, several distributions, and asymmetries are presented for 
Tevatron and the LHC at 7TeV, 8TeV, and 14TeV. 
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1 Introduction 



As the only fundamental known fermion with a mass at the electroweak scale, the 
top quark might serve as a key for understanding the fermionic mass hierarchy and is 
potentially sensitive to physics beyond the Standard Model. Since it can be produced 
copiously at the LHC in top-antitop pairs, its production cross section, its decay, and its 
properties can be studied with high precision. Evidently, accurate measurements must be 
accompanied by precise calculations. 

The first step towards precise theoretical predictions for tt production at hadron col- 
liders was made already about 20 years ago with the calculation of QCD corrections 
at next-to-leading-order (NLO) Later also electroweak radiative corrections were 

calculated [5HTU]. and recently important progress has been achieved both in the re- 
summation of logarithmically enhanced terms [TTHIl] and towards the inclusion of QCD 
corrections at next-to-next-to- leading-order (NNLO) [T51425] . First NNLO QCD results 
for the quark-ant iquark channel have been presented in Ref. [26] . 

The above-mentioned predictions are based on the approximation of stable (on-shell) 
top quarks, i.e. the top-quark decays, which proceed into pairs of W bosons and b quarks 
in the Standard Model, are ignored. A first important step towards a full NLO description 
of top-pair production and decay was made in Refs. |27l - [29] . where top-quark decays were 
treated in a spin-correlated narrow-width approximation, i.e. the top quarks are still 
on shell, but spin correlations between production and decay are taken into account. 
Meanwhile, NLO top-quark decays in narrow-width approximation were implemented in 
MCFM p]. 

Recently, first results at NLO QCD for the complete process of W"'"W~bb produc- 
tion, with intermediate off-shell top quarks and including leptonic W-boson decays have 
been obtained by two independent groups [3T1|32] . Here we give details on one of these 
calculations [31] and extend the presented results in various directions. 

The reaction pp — )• W+W~bb — )■ z/ee"'"/i~z/^bb represents one of the 2 — )■ 4 particle pro- 
cesses on the Les Houches priority list [23] • While various multi-particle NLO QCD cal- 
culations have been performed in the recent years [331 - H8] . W"''W~bb production involves 
the treatment of resonant particles for the first time in a hadron-coUider environment 
on that level of complexity. The two top resonances can be consistently treated in the 
complex-mass scheme that was introduced at the NLO level in the context of the calcula- 
tion of the electroweak corrections to the processes e"'"e~ — > W+W~ — )■ 4fermions [HUSO], 
which was the first full NLO calculation for a 2 — )■ 4 particle process. 

To compute the virtual corrections, we employ explicit diagrammatic representations 
of the one-loop amplitudes. A key feature of our approach is the factorization of colour 
structures at the level of individual diagrams. This permits to reduce the CPU cost of 
colour sums essentially to zero. Helicity-dependent structures are algebraically reduced to 
a common set of so-called standard matrix elements. In this way the sums over physical 
helicities are strongly boosted. Tensor loop integrals are related to scalar integrals by 
means of numerical algorithms that systematically avoid numerical instabilities from in- 
verse Gram determinants and other spurious singularities [^15^ . Scalar loop integrals are 
required with complex internal masses and are evaluated using the results of Refs. [53] 151]. 
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The real corrections are handled with the dipole subtraction method [551458] . and the 
phase-space integration is performed with adaptive multi-channel techniques [S5H5T| . 

Our predictions provide a complete description of hadronic W"'"W^bb production, 
including off-shell tt intermediate states, as well as contributions with a single or no 
top-quark resonance. This permits to quantify the accuracy of the narrow-top-width ap- 
proximation, which involves only on-shell tt contributions and corresponds to the Ft — )■ 
limit of our calculation, at NLO. For the integrated cross section it turns out that, in 
presence of loose cuts, the error of the narrow-top-width approximation does not exceed 
1% [2Il[S21in2]- This is perfectly consistent with the theoretically expected 0(Tt/mt) 
suppression of finite-top-width corrections to inclusive observables. 

In more exclusive measurements, such as precision rrit determinations, tt backgrounds 
to new physics that are suppressed by vetoing top resonances, or the tt background to 
H — i- lluu signals in presence of b-jet vetoes, the investigation of finite-top-width effects 
is even more important since their magnitude is not known a priori. A first systematic 
study of finite-top-width effects in exclusive observables, based on a comparison of our 
calculation against the narrow-top- width approximation of Ref. [28j, indicates that finite- 
top-width corrections to phenomenologically important observables can range from a few 
per mille to tens of per cent [02] ■ This motivates us, on the one hand, to undertake a 
more thorough comparison of the narrow- and finite-top-width approaches, which will be 
presented in a forthcoming paper [63]. On the other hand, it raises the issue of possible 
non-negligible effects resulting from the finite width of intermediate W bosons, which is 
addressed in the present paper. 

To study finite- W-width effects we consider two different descriptions of the leptonic 
W-boson decays: a (spin-correlated) narrow- W-width approximation and, alternatively, 
matrix elements for Uee'^ fi~i?fjhh production, including off-shell and non-resonant finite- 
W-width contribution^. While finite- W-width effects have already been included in the 
predictions of Ref. [32], comparing the two above-mentioned approaches provides the 
first quantitative assessment of the precision of the narrow- W-width approximation in 
i/ee'''yU~t';xbb production. In spite of the larger numerical value of Fw/iWw as compared to 
Tt/rrit, finite- W-width corrections turn out to be definitely smaller than finite-top-width 
effects. As we discuss, this is due to a double-suppression mechanism related to subtle 
cancellations between finite- W-width contributions to matrix elements and to the Ft input 
parameter. 

The paper is organized as follows. In Section [2] we discuss technical aspects of the 
calculation, including finite-top-width and finite- W-width effects (Section l2.ip . as well 
as virtual (Section 12. 2p and real (Section 12. 3p corrections. In Section 12.1.21 we also in- 
troduce a matching approach that restores exact factorization of tt production and top 
decays in the Ft — limit. This permits to avoid a few-per-cent loss of precision in 
the non-factorized NLO description of u^e'^ n~i?^hh production while including finite-top- 
width effects. In Section [3] we present numerical predictions for Tevatron and the LHC 
at different collider energies. The choice of the factorization and renormalization scales 
is discussed in some detail in Section [3^ There we introduce a dynamical scale related 
to the transverse energy of the top quarks, which avoids serious perturbative instabilities 
in the high-energy tails of differential distributions. Moreover, we advocate the use of 

^The virtual corrections are calculated in the double-pole approximation for the W resonances. 
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a collider-dependent QCD scale adapted to the particular behaviour of the leading par- 
tonic channels, i.e. quark-antiquark annihilation at the Tevatron and gluon-gluon fusion 
at the LHC. Specifically, we argue that a reduced scale provides a better description of 
W"'"W~bb production at the LHC. Results for the integrated cross sections (Section 13.31) . 
asymmetries (Section 13. 4p . and several distributions (Section 13. 5p . are supplemented by a 
numerical extrapolation of the Uee'^ fi~i?fjhh cross section to the Ft — )■ limit (Section 13. 6p . 
Finally, in Section [3?71 we show that our results are in good agreement with corresponding 
results [32] obtained with HELAC-NLO [M]. Our conclusions are presented in Section HJ 
and in App. |A]we provide benchmark results for the partonic matrix elements squared in 
lowest order and including virtual corrections. 



2 Details of the calculation 



2.1 Preliminaries 

The cross section for the hadronic process hih2 — > W^W~bb + X — > z/oe+/i~z/^bb + X 
is evaluated according to 

Cr = / dXi dX2 fa/hi{Xl) fb/h2{x2) / daab{XiPi,X2P2), (2.1) 

where the parton distribution functions (PDFs) fc/hi{xi) describe the probabilities of 
finding a parton c in hadron hi with a fraction Xi of the full hadron momentum Pi. In LO 
the colliding parton pairs are ah = gg, qq, qq, but higher orders involve contributions from 
gq, gg, gg, and gg as well. We consistently neglect contributions from initial-state bottom 
quarksjj i.e. g = u, d, c, s. Otherwise we work in the five-flavour scheme with massless 
bottom quarks. In NLO, the partonic cross section 

/ d&ab = / dao,ab + / d(Jvirt,afe + / do-^aLab + dx dafact.aft (2.2) 

J Je Je Jr Jo Je 

receives contributions from tree diagrams ((3"o,a6), virtual one- loop corrections ((3"virt,ab), 
real parton emission ((Treai.ab), and the factorization of coUinear initial-state singularities 
('^■fact.ab) into the PDFs. The subscript on the integral corresponds to the number of 
final-state particles in the reaction ab — )■ W"*"W~bb(c) — )■ z/ee"'"/i~z/^bb(c), where c is a 
possibly emitted parton. The integral over x indicates that the factorization contribution 
involves an additional integration over the fraction x by which one of the incoming parton 
momenta is reduced due to collinear initial-state radiation. 

In order to ensure the correctness of our results we have evaluated each ingredient 
twice and independently, resulting in two independent computer codes. 

In the following we discuss the treatment of top-quark and W-boson resonances as 
well as effects beyond NLO related to the factorization of tt production and top decays in 
the narrow-top- width limit. Details of the calculations of the virtual and real corrections 
are given in Sections 12.21 and | 



^We have checked that bottom-quark-induced contributions to the integrated cross section are sup- 
pressed at the level of < 0.2% at LO. 
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Figure 1: Representative tree diagrams involving two (first line), only one (second line), 
or no (last line) top-quark resonances. 

2.1.1 Treatment of unstable top quarks 

Our predictions for the process hih2 — )■ W"'"W~bb + X — i'eG'^fi~h'^hh + X provide a 
complete description of hadronic top-quark pair production and decay, including doubly- 
resonant contributions where the z/ee"'"/i~z/^bb final state results from the decay of a tt 
pair, as well as singly-resonant and non-resonant diagrams, i.e. contributions with only 
one or no top resonance. Interferences between doubly-, singly-, and non-resonant dia- 
grams are consistently taken into account. A few representative LO diagrams are depicted 
in Figured! The qq and gg partonic channels involve 14 and 31 tree diagrams, respec- 
tively, if only topologies involving two resonant W bosons are consideredjl Additional 
contributions with less than two W-boson resonances are discussed in Section 12.1.31 

To regularize intermediate top-quark resonances in a gauge-invariant way we employ 
the complex-mass scheme [50], where the top-quark width Ft is incorporated into the 
definition of the (squared) top-quark mass, 

^l = ml- imtFt. (2.3) 
Since we treat b quarks as niassless partons there are no Higgs-exchange diagrams at tree level. 
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In this way, off-shell-top contributions are consistently described by Breit-Wigner distri- 
butions, and all matrix elements are evaluated using the complex top mass /it. Technical 
implications of the complex- mass scheme at one loop are discussed in Section I2.2.6I 

The inclusive z/ee"'"/i~z/^bb cross section is dominated by the doubly- resonant top- 
pair contribution and can be described, with fairly good accuracy, in narrow-top-width 
approximation. It is thus instructive to compare our calculation to this approximation, 
which corresponds to the Ft — )■ limit. To avoid confusion between the treatment of top- 
quark and W-boson decays, in the following we denote the Ft — and Fw — ^ limits as 
narrow-top-width (NtWA) and narrow- W-width (NwWA) approximations, respectively. 
Contributions that vanish in NtWA and NwWA are called finite-top-width (FtW) and 
finite- W-width (FwW) effects, respectively. Our treatment of FwW effects is discussed in 
Section I2.1.3[ 

For what concerns top resonances, we point out that FtW contributions are included 
everywhere in our calculation, i.e. we never make use of the NtWA. Nevertheless, in the 
following we briefly introduce this approximation in order to discuss the origin of FtW 
effects and other features of our predictions. In the NtWA, each top-quark resonance 
leads to 

1 vr 

1™ 0^ oT^ = — ^ ^(Pt ~ "^t)? (2-4) 

where the 5-function that forces the top quark on its mass shell is accompanied by a 
1/Ft factor. In NtWA the z/ee"'"/i~z/^bb cross section includes only contributions involving 
two top resonances, which are proportional to 1/F^. Singly- and non-resonant diagrams, 
as well as their interference with doubly-resonant diagrams, are neglected due to their 
suppression in the Ft — ?■ limit. As a result of these approximations, the differential 
hih2 — 7- tt — )■ z/ee^/i~t/^bb cross section is factorized into the hih2 — )■ tt production cross 
section times the t — )■ Wb — )■ /z//b partial decay widths, 

d(TNtwA = Fj"^ (do-tt dFt_,^ dF^_^j.) , (2.5) 

where the subscripts i,j refer to the (anti)top-decay final states u^e'^h and fi^u^h, and 
the total top-quark width is obtained by summing over all relevant decay channels. 



dr.^i. (2.6) 



Top-quark spin correlations in fl2.5p are implicitly understood. 

In NtWA, LO and NLO partonic cross sections can be schematically written as 

d<wA = (r[^^)-^[d-t^°ciF,^OdF,^«. 

da?Jt"^^A = (rr°)-^ [(d4 + d4) dF^, dF^, + d^ (dFt^, dF,%. + dF°^, dF,^,.) 

(2.7) 

where the superscripts and 1 indicate LO and correction contributions to NLO quan- 
titles, i.e. dFr_^o ^ ^F^, + dF^^^, and da^^^o ^ ^^o + ^ote that d^ ^ da^^^, 

since the ingredients of dcrJj^A ^^d dcrJJ^wA have to be evaluated with input parameters 
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at the corresponding perturbative order, i.e. da^f and dcr^^ are evaluated with PDFs and 
as in LO and NLO, respectively. Since the top decay does not involve as at LO we have 
dr°_j.j;. = dr}^°^. Note also that LO and NLO predictions must be computed using total 
decay widths Fji"^ and rf'"^, respectively. Naturally, this also holds for our calculation, 
not only in NtWA. 

As a result of the truncation of the perturbative expansion, the NtWA of the NLO 
cross section (12.71) involves only three contributions, where the corrections are applied 
either to dati or to one of the decays. A pure fixed-order NLO calculation does not 
include contributions like dal^dTl^^dT^^j, where the NLO correction is simultaneously 
applied to the production and decay parts of the process. Such contributions are formally 
of NNLO, but given their non- negligible numerical impact it is desirable to take them 
into account. As discussed in Section [2.1.21 this can be achieved by means of a simple 
prescription which is applicable also in presence of FtW and FwW corrections and is 
derived by matching Eq. fl2.7p to Eq. fl2.5p at the level of the fully inclusive cross section. 

As compared to the NtWA, our calculation includes two types of additional effects ow- 
ing to the FtW: corrections resulting from the top-quark off-shellness in doubly-resonant 
channels and contributions from singly- as well as non-resonant diagrams. At NLO, FtW 
effects manifest themselves also in the form of non-factorizable corrections to doubly- 
resonant diagrams, which originate from one-loop topologies or interferences between 
real-emission diagrams, where top-antitop production and top decays are linked via ex- 
change of QCD partons. Technically, the calculation of non-factorizable, singly-, and 
non-resonant one-loop contributions involves pentagon and hexagon diagrams, which rep- 
resent a much higher level of complexity as compared to the four- and lower-point diagrams 
that appear in NtWA. 

The non-factorizable contributions of virtual and real origin are enhanced by large 
logarithms of Ft/mt originating from soft gluons. However, it is well known that — in suf- 
ficiently inclusive tt observables — such logarithms cancel, and the remaining FtW effects 
yield rather small contributions of order Ft/m^ ^ 0.9% [53H^ . Nevertheless, a reliable 
quantitative assessment of FtW effects is important in order to achieve per-cent-level 
precision in the (inclusive and differential) description of tt production. Comparing our 
predictions to the NtWA, one can obtain a precise determination of FtW contributions to 
any infrared-safe observable. For the case of the integrated cross section, we performed 
this comparison by means of a numerical extrapolation of our predictions to the Ft — )■ 
limit [21] • The results, discussed in Section [XBl demonstrate that the NtWA for the 
integrated cross section agrees with the full calculation at the sub-per-cent level. 

As already mentioned in the introduction, in more exclusive phenomenological studies, 
such as precision rrit measurements or Higgs- and new-physics searches, the size of FtW 
corrections cannot be anticipated a priori and turns out to range from a few per mille to 
tens of per cent [S21E3]- 
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2.1.2 Matching to NLO inclusive tt cross section 

Let us now discuss effects related to the truncation of the perturbative expansion at 
NLO in the presence of unstable intermediate particles. To start with, we consider the 
fully inclusive cross section in NtWA, 

dfTNtwA = o"tt BR^_^j BRj^^-, (2.8) 

which is obtained by integrating ( 12. 5 p over the full phase space and is given by the on-shell 
inclusive tt cross section, 



cTtt = J d(Ttt, (2.9) 

times the branching fractions 

BRt^. = = (2.10) 

i t J- 1 

with k = i,j. Apart from Coulomb effects near threshold, the above relation between 
the pp — J- tt — )■ ij and the on-shell tt cross sections is valid to all orders of perturbation 
theory However, due to missing higher-order terms, the NLO approximation (12.71) 
does not fulfil (12.81) exactly. The mismatch can be expressed in terms of products of NLO 
contributions as follows, 



A^NLO NLO -RoNLO -RoNLO /^^NLO 

= (rf^°)-^ \K + 4) rU. + 4 (rU r^,, + r^, r^^,)! . (2.11) 



^NLO 
'NtWA 



Rewriting (12. lip as a relative correction to (12. 8p yields 

A NLO 

s^NLO _ ^'-'trunc 

a,fOBRf_i?BRfi;,^ 



^trunc - .NLORpNLORpNLO ~ [^«(^ + ^tt + ^iXj, (2.12) 



where the factors 



and 



Stt = 


4 


= 1 - 


4 




4^° 


Xk = 


t-s.fc 
pNLO 


= 1 - 


T-O 

Wfc 
pNLO 



(2.13) 



(2.14) 



represent NLO corrections to tt production and decay, respectively. For the case of a 
di-lepton final state, where Eq. (I^TT^ simplifies to 



Cnc = 2x(l-x)5tt + x^ (2.15) 

Since x ~ 10% and 6tt — 10—30%, the correction 5jJu^ can reach a few per cent and 
should thus be taken into account. 
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In the case of di-leptonic decays of the tt system, this can be achieved by using the 
approach of Refs. ISSUSO], where the factor (rNLO)-2 ^ (jpo ^ ^1^-2 (j^^j) is expanded 
and truncated at NLO. The corresponding differential NLO cross section reads 



d^e^r = (r?)"^ (d4 + d4)dr^,dr 







da°- ( dr^ dr^ + dr° dr-^ - 2 — dr° dr^ 



(2.16) 



Integrating over the full phase space, and expressing dr°_j.^ and dr^_j.^ in terms of the 
total top width and W-decay branching fractions via 



t po ■ 
J- w 



pi pi 

-L + -L T 



W->fc 



pi 

J- w 



(2.17) 



one obtains the following expression, where the total and partial top-decay widths cancel 
out, 



da. 



NLO 

exp 



4 + 4 E 

fc=j,J 



T} 



po 



W-»fc 
pO 



pO 
J- w 



n 



pO 



n 



pO 

po • 

J- w 



(2.18) 



Note that, consistently with the treatment of the top width in (12.161) . we expanded the 
NLO W-width term in (l2Tni as = - 1^/(1^)2. Moreover, to be fully 

general, in (I2.17P we also included corrections r^_^^ to the W decays. Comparing fl2.18p 
to the factorized expression (12. 8p with NLO branching fractions 



BR 



NLO 
't->fc 



pNLO 
-*- t-^k 

pNLO 



Ff^^BR 

pNLO 



NLO 
W->fc 



pO 
J- w- 



W-j-fc 



pO 
J- w 



w 



(2.19) 



yields 



NLO 

tt 



750 NLO -n-QL^ 



)NLO 



j^NLO ^1 



pO 

n " 



W->fc 



pO 
J- w 



W 



pO 

n " 



W-j-fc 



pO 



+ r", 



J- w- 



pO 
J- w- 



w 



pO 
J- w 



pO 

n " 



pO 
J- w 



n 



pO 



pO 



(2.20) 



This indicates that, in general, the inclusive cross section resulting from (I2.16P is not 
identical to (12. 8p . The mismatch (I2.20p is formally of NNLO and is due to the corrections 
Tw^fc to the exclusive W decays and to the fixed-order expansion of l/(r^ + L^). Thus, 
if the W decays do not receive NLO corrections the fixed-order approach (I2.16P can be 
reconciled with (12. 8p by avoiding the fixed-order expansion of the W-width term, i.e. using 
l/^w^^ = l/(rw + rw) everywhere. This corresponds to replacing l/F^ ^ i/(ro^+riv) 
and omitting the Ty^/T'^ term in (12.201) . 
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In practice, this is applicable to di-lepton final states, but not to hadronically decaying 
top quarks or in presence of electroweak corrections. Therefore we prefer to ensure the 
validity of f l2.8p at NLO in a different way. Instead of removing higher-order contribu- 
tions by expanding and truncating the term (Ff which originates from the Dyson 
resummation of the top-quark self-energy in (12. 4p . we keep the exact NLO top-quark 
width everywhere and supplement our calculation by the missing higher-order contribu- 
tions (12. lip . This is done by correcting the normalization of the NLO pp — )■ tt — )■ ij cross 
section in NtWA, 

d«: = (l + 5^^°")d<o^, (2.21) 

with a matching factor 

xNLO 

rNLO+ _ "trunc fn nn\ 

" ~ 1 ^NLO yz.zz) 

^trunc 

that restores the consistency with the on-shell tt total cross section, i.e. 

For the total cross section, the approximations (I2.16P and (I2.2ip are fully equivalent (for 
di-lepton final states and if F^^^ is used everywhere). However, they yield different 
predictions for exclusive observables, in which case we expect (I2.2ip to be more accurate, 
thanks to the presence of higher-order contributions. 

The matching correction (I2.22p guarantees a consistent inclusive cross section also in 
presence of FtW (and FwW) effects. In this case, for the corrected cross section 

da'^'^T^- = (1 + 5^^°^) da^L? _- ,r, (2.24) 

we have 

which remains consistent with (12. 8 p for Ft — 0. Moreover, the matching procedure does 
not disturb the 0(Tt/mt) finite-width contributions Acxp^^ to c"i,^e+/i-i/^bb) since the inter- 
play of these two corrections represents a suppressed NNLO effect of order 5'^'"*-*^Ft/mt. 

Finally, let us note that in order to determine the correction factor (I2.13P that enters 
^NLO+^ instead of performing an explicit on-shell tt calculation we use 



< /d^^cuts ^ /d<eV-P.bb^cuts /F 



■^NLO 
t 



(2.26) 



i.e. we replace the fully inclusive a^i cross section by the u^e'^ fi~u^hh cross section in 
presence of the (fairly inclusive) cuts specified in Sectional By definition, the various 
(T° and a^^^ cross sections in f l2.26p must be evaluated using NLO PDFs, NLO a^, and 
F[^'"'-', and the factor (Ff^°/F}i"°)^ on the r.h.s. of (I2.26P compensates for the fact that, 
in order to match the l.h.s, d(j° g+^_p should be computed with a LO top-quark width. 
We expect the above approximation to be quite precise since the mismatch between the 
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tt and z/ge"''/! u^hh NLO cross sections in (I2.26P induces a strongly suppressed error of 
order (5Nlo+)2 

in a^^^ __ , r. Moreover, cut and finite- width effects in fl2.26p cancel to a 
large extent in the ratios and are further suppressed by the factor 2a;(l — x) in (12.151) . 

In the derivation of the above NLO'*' correction factor we implicitly assumed that the 
top-quark width input parameter and the matrix elements are evaluated at the same QCD 
scale. However, while rrit represents the natural scale choice to compute F^^*-*, in general 
the matrix elements for i/^e'^ fi~Ufjhh production might be evaluated using a different scale 
/X, like the dynamical scale proposed in Section 13. 2[ which adapts to the hard scale of the 
tt production part of the process. Since the scale fi also enters the top-decay part of the 
matrix elements, a scale choice fi ^ rrit implies that partial and total top-decay widths 
are evaluated at different scales. This mismatch can be compensated by redefining the 
NLO differental cross section as 



da^LO 



rr-M Y ^^NLO (2.27) 



PWC V rf^°(/") / ^-'=+/^-^Mbb' 



where the overall factor {r^^'~'{mt) /T^^'^lfi))'^ effectively restores /i — )■ mt in the partial 
decay widths. The (dynamical or fixed) scale n is the one used in the matrix elements, 
including variations in scale-dependence studies. The relative shift of the NLO cross sec- 
tion induced by the correction (I2.27P amounts to about —0.02 x (/i/mt — 1) for fi/rrit ~ 1. 
In principle this correction can be applied to all NLO predictions. However, we decided 
to present fixed-order NLO results in a more conventional way (scale /x in matrix elements 
and rf^'^ at the fixed scale rrit) and to include the above correction only in NLO^ pre- 
dictions. To this end, we evaluate fl2.24p and fl2.26p using "partial-width corrected" NLO 
predictions (I2.27p . In presence of a dynamical scale, /i = /Xdyn, instead of a phase-space- 
dependent correction we apply a constant factor Tf^'^{mt)/Tf^^{fidyn) evaluated at the 
logarithmic average of /^dyn (see Table [1]) . 

2.1.3 Treatment of unstable W bosons 

To describe the leptonic W-boson decays, — )■ u^^^ and — t- ii^u^, we employ 
two different approaches: a spin-correlated NwWA and, alternatively, the full set of di- 
agrams contributing to i/ee'''/i~t'^bb production. Analogously to the case of top-quark 
decays, the NwWA includes only contributions with two on-shell W bosons, while the 
FwW calculation involves additional effects from off-shell W bosons and from singly- 
resonant u^e'^n'i'^ production]^ In NwWA, the effect of the W-boson width is retained 
only in the W propagators — approximated by vr/ (Mwr-w)^(Pw ~ ^w) — while elsewhere 
we set Fw = Fz = and use real- valued parameters Mw and Mz- In the alternative calcu- 
lation with the full set of diagrams, to include FwW effects we employ the complex-mass 
scheme, i.e. we use complex W- and Z-boson masses and a corresponding complex-valued 
weak mixing angle, 

/i^ = - iMwFw, /i| = M|-iMzFz, 0^ = 0^^=—. (2.28) 



*A11 channels contributing to v,,^'^ i^/ibb production involve at least one resonant W boson. 
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Figure 2: Examples of LO diagrams that involve a single W-boson resonance and con- 
tribute to the FwW corrections. 

Examples of tree diagrams in NwWA and additional singly-resonant diagrams con- 
tributing to the FwW calculation are shown in Figure [T] and Figure Ej respectively. In 
all singly-resonant diagrams, the four-lepton final state originates from Z,7 — )■ z/ee"'"yU~z/^ 
subtopologies involving one intermediate s-channel W"*" or W~ propagator. We note that 
such topologies do not involve any top-quark resonance and are thus expected to be 
strongly suppressed. Their inclusion increases the number of LO diagrams in the qq (gg) 
channel from 14 (31) to 38 (79). 

At NLO, a resonant W-boson pair decaying into z/ee"'"/i~z/^ does not receive (factoriz- 
able or non-factorizable) QCD corrections. This renders the inclusion of FwW effects fairly 
straightforward, as compared to the case of FtW effects. Since FwW contributions are 
strongly suppressed already at LO (see Section [3]), we adopt a double-pole approximation 
to include FwW effects in the infrared-finite part of the virtual QCD corrections, which 
is defined by the sum of the virtual corrections and the contribution of the /-operator 
of the real corrections in the dipole subtraction approach [551158] . This approach, which 
is discussed in Section I2.2.7[ involves only doubly- W-resonant diagrams. Similarly as in 
NwWA, the matrix elements in double-pole approximation are evaluated with on-shell 
W bosons and using Fw = Fz = and real W- and Z-boson masses everywhere, apart 
from the W propagator. The latter reads l/(pw ~ + ifw^w) and takes into account 
off-shell-W effects. Apart from the (finite part of the) virtual corrections, in all other NLO 
contributions FwW effects are included exactly, using the complex- mass scheme (12. 28 p . 

For what concerns the expected magnitude of FwW contributions to inclusive tt ob- 
servables, one might naively expect an impact of order Fw/A/w — 2.5%, i.e. three times 
larger than the typical size of FtW effects. However, more careful considerations indicate 
that the consistent inclusion of FwW corrections in matrix elements and input parame- 
ters leads to doubly-suppressed contributions of 0{ ^^^^ ) in the inclusive cross section. 
This is due to the fact that, in the Ft — )■ limit, the only contribution to the inclusive 
VeO^ li~Vyhh cross section (I2.25p that involves W-boson decays, and is thus sensitive to 



Fw, is the branching ratio 

-OpNLO _ J "^-^ t-¥lvih /r, r)Q\ 



Moreover, owing to the identical Fw dependence of the numerator and denominator of 
( I2.29P — which correspond to the integrated top-decay matrix elements and the Ft input 
parameter in our calculation — FwW corrections cancel in the branching ratio. Thus the 
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consistent inclusion of FwW corrections — in the ah — ?■ z/ee^/i^z/^bb(c) matrix elements 
and in the calculation of the Ft input parameter — does not affect the leading contribution 
to fl2.25p and yields only tiny corrections of to the fully inclusive cross section. 

These considerations do not depend on the inclusion of the matching correction (12.241) 
and remain valid also in presence of contributions of type (12. lip , which violate (12. 8p . We 
also point out that omitting FwW effects in the top-quark width used as input parameter 
would induce a fake FwW shift of 0(Fw/Mw) in the integrated cross section. 

The validity of the double-suppression mechanism is restricted to the case where the 
top-decay phase spaces are fully integrated over, and the presence of cuts is expected to 
lead to additional FwW effects resulting from the incomplete cancellation of (9(rw/Afw) 
corrections to the numerator and denominator of (I2.29p . However, as long as the cuts are 
rather inclusive, also such FwW effects are expected to remain well below C(rw/Mw). 
This justifies the approach of Ref. [21], where off-shell and non- resonant contributions of 
the top quark to W"'"W~bb production were combined with W-boson decays in NwWA. 

For very exclusive observables, in contrast, it is important to investigate whether 
FwW corrections can become non-negligible. To this end, in this paper we study FwW 
effects by comparing several distributions obtained in NwWA and with the FwW variant 
of our calculation. While FwW effects have already been included in the predictions 
of Ref. [22], our study provides the first quantitative assessment of the precision of the 
NwWA in i/ge^/W ^'/xbb production. 

2.2 Virtual corrections 

2.2.1 Diagram-by- diagram approach 

The sum of the LO and virtual contributions to the partonic cross section is derived 
from the LO and one- loop matrix elements Aio and Aii according to 



jdao,ab + jda,M,ab = ^ j d<l>6 | J] ^ |Mo,a6p + XI ^ ^^(A^l.ab -Mo.ab) }> 
^ ® col pol col pol 

(2.30) 

where the squared partonic centre-of-mass (CM) energy s appears in the overall flux factor 
1 / (2s) and d$6 is the phase-space element of the six-particle final state. The sums run 
over all colour and helicity states and implicitly include the averaging over initial-state 
colours and helicities. We calculate the amplitudes M.q and M.i for qq/gg — )■ W'^W~bb — )■ 
i^ee'^ IJ,~u fjhh in terms of explicit Feynman diagrams and algebraically reduce the one-loop 
diagrams to spin structures ("standard matrix elements") and standard one-loop tensor 
integrals, which are subsequently evaluated numerically, as described in more detail below. 

The qq and gg channels comprise 294 and 795 one-loop diagrams, respectively (see 
examples in Figure The most complicated ones are the 96 pentagons and 21 hexagons 
that contribute to the gg channel and involve tensor integrals up to rank five. Feynman 
diagrams are generated with two independent versions of FeynArts [681169] , and one- 
loop amplitudes are reduced as already described for pp — )■ ttbb in Refs. [70l[7T] using two 



While our calculation includes the contributions of all three quark generations, we count diagrams 
with closed fermion loops only for a single generation. 
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Figure 3: Examples of one-loop diagrams contributing to qq /gg — W'^W^bb — ?■ 
z/ee"'"/i~z/^bb: doubly-top-resonant diagrams with corrections to tt production or decay 
(first line), non-factorizable pentagons and hexagons with two top-quark resonances (sec- 
ond line), pentagons and hexagons with less than two top resonances (third line). 

in-house Mathematica programs, one of which relies on FormCalc [72] for preliminary 
manipulations. 

The employed approach strongly mitigates the complexity inherent in Feynman dia- 
grams by exploiting factorization of colour matrices, reduction of helicity structures to 
compact spinor chains, and recycling a multitude of common subexpressions. The re- 
duced expressions are automatically converted into Fortran?? programs that evaluate 
colour/helicity summed quantities with very high CPU efficiency. 

The virtual corrections are obtained from the interference of the one-loop and LO ma- 
trix elements summed over external-state colours and helicities on a diagram-by-diagram 
basis. 

2. 2. 2 Colour factorization 

One of the key features of the diagram-by- diagram approach is that the cost related 
to the large number of diagrams is compensated by the possibility to perform colour sums 
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very efficiently. This is a consequence of colour factorization: individual (sub) diagrams, 
r, consist of a single colour-stripped amplitude A^^^ multiplied by a single colour factor 

=^(r)c(r). (2.31) 

More precisely, each diagram gives rise to 3"* colour-factorized contributions of type 
fl2.3ip . where is the number of quartic gluon vertices in the diagram. These terms 
are handled as separate subdiagrams. However, most diagrams do not involve quartic 
couplings, and their colour structures factorize completely. The colour factor C^^^ can be 
linearly decomposed in terms of a fixed basis of colour structures Ck, 



C^'' = j:crC„ (2.32) 



k 

(T) 

where the coefficients c), are just a set of numbers for each diagram F. For the qq and 
gg channels there are only 2 and 3 different colour structures, respectively, which may be 
taken to be 

Cf- = ^ab^c^c^, Cf- = ifabcT^^c^, Cf^ = dabcT^^c^. (2.33) 

Here a, b and Cj denote the colour indices of the two incoming gluons and the fermion i, 
respectively, are the group generators in the fundamental (quark) representation, and 
fabc, dabc s^TQ the usual SU(3) constants. 

The summation over external colours is performed once and for all at the level of the 
colour basis and the LO matrix element. To this end, we compute the colour-interference 
matrix 



Y^CkCl (2.34) 



col 

and, reducing the full tree matrix element in colour space, 

Mo = $^A^f°^C;, (2.35) 

we build the interference of A^o with the elements of the colour basis as 

Mf"^^ = Y.C,Ml = Y^h\Mr^y . (2.36) 

col I 

Then, upon reduction of the factorized colour structure of the loop diagrams, 

^(r) ^ ^(r)c(r) ^ ^(r) ^ip^^ ^ ^2.37) 

we obtain the colour-summed interference between A^q and the complete one- loop ampli- 
tude Mi = Y.r^^^^ as 

Y^M.Ml = J2J2^^^'^l = J2A^^^ (j:cPMr'] ■ (2.38) 

col col r r \ A; / 
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The colour-summed result is given by a combination of previously computed colour-Born 
interference terms fl2.36p . For each phase-space point, this requires a single evaluation of 
the non-trivial colour-stripped amplitude A^^'^ of each (sub) diagram. 



2.2.3 Algebraic reduction of helicity structures and helicity sums 

The helicity structures encountered in the explicit evaluation of all Feynman diagrams 
are algebraically reduced to a common basis of standard matrix elements (SMEs). The 
general form of SMEs for the a{ki)h{k2) — ?■ W'^(A;3)W~(fc4)b(fc5)b(fc6) channel for the 
initial states ah = qq /gg is 



X 



M^lr = Q::i:^r:...4^r(^i)<'(^2)5^:(fc3)er (^4) [%{k,)i., . ..^..uru^ik,)] , (2.39) 

where Q^-'fisl^ Qm-'iliLfj.4 consist of combinations of metric tensors and external mo- 
menta, and a, r = ± refer to the chirality projectors u± = (1 ± 75)/2. In the double-pole 
approximation (see Section l2.2.7p . W-boson decays are described via effective polarisation 
vectors 



^w+(^3) 



eu{ky^)^^uj.v{kc+] 



\/2s„((fc,^ + k,+ f -M^ + iMwF 



w 



V2s„(^(V + -M^ + iMwFwj 

which include the left-handed lepton currents and the W-boson propagators. In our cal- 
culation we encounter about 800 and 2000 SMEs for the qq and gg channels, respectively. 
These compact spinor chains permit to decouple helicity information from the remnant 
parts of the diagrams, so that helicity sums can be performed in a diagram-independent 
and efficient way. In practice, the colour-stripped part of each loop diagram [see (I2.37P 
is expressed as a linear combination of SMEs and tensor integrals, 

m 

N-l 

J"^^ = E E ^'Sn..,.^;!..,, + rational parts, (2.41) 

P ju-Jp=0 

(T) 

where the index m here also includes the chirality indices a, r. The coefficients IC^'j^ jp are 
rational functions of the kinematic invariants. These functions involve only denominators 
from intermediate-particle propagators and are free of spurious poles that might generate 
numerical instabilities. The functions Tj^ are the coefficients of the one-loop tensor 
integrals, whose evaluation is briefly described below. 
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Helicity sums are performed at the level of the interference of the diagram-independent 
SMEs with the colour-projected Born amplitude f l2.36p . 



= Y^MmMf""^ = J2lkiJ2^m (2.42) 

hel I hel 

This matrix is computed only once per phase-space point employing the Weyl-van der 
Waerden spinor formalism of Ref . [73] . Using Mkm one can directly obtain the colour- and 
helicity-summed contributions of each loop diagram in terms of its colour- and helicity- 
independent form factors J-m and the coefficients of its factorized colour structure 

can), 

col hel col hel r F m \ k / 

2.2.4 Covariant decomposition and numerical reduction of tensor integrals 

Tensor one-loop integrals with N propagators and P Lorentz indices are expressed in 
terms of totally symmetric covariant structures {g . . . gp . . ■p}jl "j'p involving g'^'^ and the 
external momenta Pi, ■ ■ ■ ,Pn-i, 

(27r/i)^-^ f ^ q^-...q^^ ^ ^ 

1^ J 11^=0 [(g + P»)2-< + iO] j,„^=o 

(2.44) 

with D denoting the number of space-time dimensions. For details of the notation we refer 
to Ref. [52]. To describe iV-point integrals with > 5, tensor structures with only four 
external momenta would be sufficient. However, in order to avoid potential instabilities 
due to inverse Gram determinants we use a redundant set of structures, including the 
metric tensor and — 1 momenta. 

The virtual corrections to qq — )■ W"'"W~bb and gg — )■ W"'"W~bb involve tensor integrals 
up to ranks P = 4 and P = 5, respectively. As sketched above, the one-loop amplitudes 
are expressed as linear combinations of tensor-integral coefficients T^^ The latter are 
evaluated by numerical libraries that recursively reduce them to master integrals using 
the methods of Refs. [5ll|52]. Avoiding an exphcit reduction of analytic expressions to 
master integrals, this numerical approach prevents prohibitively large expressions and 
permits to adapt the reduction strategy to the specific numerical problems that appear 
in different phase-space regions. 

Tensor A^-point integrals with A^ > 5 are expressed in terms of lower-rank and lower- 
point integrals exploiting the four- dimensionality of space-time [5T|[52] J§ The tensor rank 
and the number of propagators are simultaneously reduced without introducing inverse 
Gram determinants. Consequently, the maximal power of inverse Gram determinants 
resulting from the entire reduction is given by the maximal rank of four-point integrals, 
which never exceeds four in renormalizable gauges. Scalar hexagons and pentagons are 
reduced to boxes using Melrose's method Tensor 4-point and 3-point integrals are 



'Similar reductions are described in Refs. [7il[75] . 
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reduced to scalar integrals with the Passarino-Veltman algorithm [77] as long as no small 
Gram determinant appears in the reduction. If small Gram determinants occur, alterna- 
tive schemes are applied More precisely, we make use of expansions of the tensor 
coefficients about the limit of vanishing Gram determinants and possibly other kinemat- 
ical determinants. One- and two-point tensor integrals are obtained with numerically 
stable analytic expressions. 

Ultraviolet (UV) divergences are regularized dimensionally throughout, but infrared 
(IR) divergences are treated in different variants, which comprise pure dimensional regu- 
larization with strictly massless light quarks and a hybrid scheme with small quark masses 
and massless gluons. The corresponding scalar integrals are evaluated using the methods 
and results of Refs. [53l|5ll[80], and different regularization schemes are translated into 
each other as described in Ref. [81] . 

The calculation of tensor integrals is implemented in two independent Fortran li- 
braries. This permits to perform detailed cross checks, which confirm the excellent nu- 
merical stability of the reduction procedure. An automatic cache system is implemented 
that strongly boosts the reduction by recycling a multitude of tensor integrals among 
Feynman diagrams with common sub-topologies. 



2.2.5 Rational parts 

In D = 4 — 2e dimensions, UV-singular tensor integrals give rise to 1/euv poles and 
UV- finite remainders Tj^,,,jp, 

Consequently, their D-dependent prefactors /(-D) need to be expanded in D — 4, 

= - 2f (4X..,,, (2.46) 

resulting in so-called rational terms that are proportional to the pole residues R^^ j^. 
Rational contributions originate from D-dependent terms in tensor-reduction identities 
and in the loop-momentum-independent part of the diagram numerators. 

We employ the treatment of rational terms of ultraviolet or infrared origin as described 
in Appendix A of Ref. [70]. We use the fact (proven in Appendix A of Ref. [70]) that in 
the 't Hooft-Feynman gauge and similar gauge fixings rational terms of IR origin cancel 
in truncated one-loop amplitudes and only contribute via wave-function renormalization 
factors. Rational terms of UV origin are obtained in a straightforward way by performing 
the relevant expansions about D = 4 automatically by means of a catalogue of residues 
of UV poles (see e.g. Ref. [5^). 



2.2.6 Treatment of unstable top quarks 

The presence of intermediate unstable top quarks in hih2 — )■ W"'"W~bb — )■ Uee'^ fi'u^hh 
represents a non-trivial new aspect as compared to previous NLO QCD studies of multi- 
particle processes. To regularize intermediate top-quark resonances in a gauge-invariant 

^Similar procedures based on numerical evaluations of specific one-loop integrals |74p78j or expansions 
in small determinants ^79j have also been proposed by other authors. 
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way, we employ the complex- mass scheme [50j. In this approach the top-quark width Ft 
is incorporated into the definition of the renormalized (squared) top-quark mass fl2.3p . 
This complex parameter /i^ is identified with the position of the pole of the top-quark 
propagator, and the top-mass counterterm 5/it is related to the top-quark self-energy at 



Pt = l^t via [50] 



SfM, 



[S*'V?) + S^.L(^?) + 2S*'^(/.,^)] 



(2.47) 



where S*'^, S*''", and E*'^ are the left-handed, right-handed, and scalar parts of the top 
self-energy, respectively. This yields the one-loop counterterm 



TT 



+ 



4 



(2.48) 



in Z) = 4 — 2e dimensions, where n is the reference scale of dimensional regularization. The 
evaluation of one-loop scalar box integrals in the presence of complex masses represents 
another non-trivial aspect of the complex-mass scheme. In our calculation we employ 
the results of Ref . [M] , where explicit analytic continuations of box integrals have been 
presented for all kinematic configurations that are relevant for physical processes. 



2.2.7 On-shell projection for off-shell W bosons 

As discussed in Section 12.1.31 we treat the leptonic W-boson decays in (spin-correlated) 
NwWA and, alternatively, including also FwW effects. In the latter case, we employ 
exact LO and real-emission amplitudes — including all non-resonant and off-shell effects 
resulting from the FwW — together with the one-loop matrix elements and the /-operator 
from dipole subtraction in "double-pole approximation" (DPA). The DPA represents the 
leading contribution in an expansion of the corresponding matrix elements around the 
resonance poles in the W-boson propagators. 

The pole expansion is based on the idea [82l[83] to separate the resonant part of an 
amplitude Ai which peaks at ~ in the following way, 

= ]^ + ) = + ,Lm. + ^2.49) 

where R{M'^) is the gauge-independent residue of the resonance and the non-resonant 
contributions from R{p'^) — R{M'^) and N{p'^) are neglected in DPA. A consistent im- 
plementation of this idea, however, involves some complications, in particular, if two 
resonances instead of one are involved. Nevertheless this concept was very successfully 
applied, for instance, to W-pair production in e"'"e~ annihilation at LEP2, as e.g. reviewed 
in Ref. ^1]. Here we follow the "hybrid concept" of Ref. [85j, where only the IR-finite 
virtual corrections are treated in DPA, while keeping the lowest-order and the remaining 
real-emission matrix elements fully off shell. One of the mentioned subtleties concerns 
the appearance of so-called non-factorizable corrections, which involve non-analytic terms 
like ln(p^ — M^) in the resonance region. Such terms are ignored in the naive equation 
fl2.49p : they are caused by soft photon or gluon exchange between production and decay 
processes that are linked by the resonance. In our case, such non-factorizable corrections 
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do not appear, since we consider QCD corrections in combination with purely leptonic 
W-boson decays, i.e. there is no gluon exchange between W production and decay. Thus, 
only so-called factorizable corrections are relevant for us, which are just the corrections 
to the residue R{M'^). 

In this context, a second subtlety arises in the evaluation of R{M'^) / {p^ - M^) concern- 
ing the kinematics. The residue R{p^) typically depends on the momenta of all involved 
particles, and the consistent substitution — )■ requires a deformation of momenta 
on the full off-shell phase space of the reaction. This "on-shell projection" is needed in 
order to define a gauge-invariant set of corrections and maps each phase-space point into 
an associated phase-space point with on-shell W bosons. This procedure actually involves 
some freedom, but this ambiguity changes the results only at the level of non-resonant 
terms. 

Technically, we proceed as follows to obtain the virtual corrections in DPA for the 
W-boson resonances in qq/ gg — ?■ z/ee'^yU~z/^bb. We select those diagrams that involve 
two resonant W bosons and perform a projection of the final-state momenta that puts 
W bosons onto their mass shells. Since we deal with corrections of 0{as) in DPA only, the 
intrinsic ambiguity leads to differences of the order of asT^f^/{'I^M^J^) for different versions 
of the projection, provided that the projection does not induce large changes in the matrix 
elements. To this end, we must keep the invariant masses of the top and antitop quarks 
fixed. In fact, we do not modify the four-momenta of the top and antitop and define the 
on-shell projected momenta pi in terms of the original momenta pi as follows: 

Pw+ =Pt- Ph, 

Pv, =PW+ -Pe+, 

Pw- =Pt-Ph, 
Ppm=Pw--Pm-- (2-50) 

Each top decay is treated separately, the momenta of the bottom quarks and leptons are 
rescaled such that = and Pl^ = = p?^, and the momenta of the neutrinos 
and W bosons are determined from four-momentum conservation. The above-mentioned 
ambiguity in the on-shell projection, e.g., manifests itself in our choice to rescale the 
charged-lepton momenta and to fix the neutrino momenta via momentum conservation; 
the analogous on-shell projection with the roles of charged leptons and neutrinos inter- 
changed would be equally good. Apart from the W propagators, the matrix elements are 
computed using on-shell projected momenta, physical (i.e. real) W- and Z-boson masses, 
and F-w = Tz = 0. The phase space and the W propagators are evaluated using the 
original off-shell momenta and including the W-boson width. 

This procedure accounts for the most important off-shell effects in spite of the use of 
the DPA for the one- loop matrix element. This fact is strongly supported by the results 
on W-pair production in e"'"e~ annihilation, where comparisons between DPA and full 
off-shell calculations exist for LO [HS] and NLO electroweak corrections |^I50]. 



Ph = Ph^ , 

Pe+ = Pe 



2pw+Pe+ 



Ph = Ph- o 
Pi,- = Pi, 



'2pW'Pa- 
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Figure 4: Examples of real-emission diagrams with two (first two lines), one (third line) 
or no (last line) top-quark resonances. All depicted diagrams involve two W-boson reso- 
nances. 

2.3 Real corrections 

2.3.1 Matrix- element evaluation 

The real corrections receive contributions from the 2 — )■ 7 partonic processes gg — )■ 
z/ee"'"/i~z/^bbg, qq — )■ v^^'^ n'u^hg, gq — )■ u^e'^ fi'u^hhq, and gg — > u^e'^ ^~u^hhq. In 
the NwWA the gg channel involves 208 tree diagrams, while the qq, gq, and gg channels, 
which are related by crossing symmetry, are described by 90 tree diagrams each. Examples 
of doubly- W-resonant real-emission diagrams are depicted in Figure |H When including 
FwW effects, additional singly- W-resonant diagrams must be taken into account (for 
examples see Figure [5]) resulting in a total number of 508 and 234 diagrams for the gg 
channel and the qq channel, respectively. The 2—7-7 matrix elements are evaluated 
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Figure 5: Examples of real-emission diagrams that involve a single W-boson resonance 
and contribute to the FwW corrections. 



with Madgraph |87] and, alternatively, using the Weyl-van-der-Waerden formalism of 
Ref. [73]. In addition, in presence of FwW effects we used OpenLoops [88] to evaluate 
the real-emission matrix elements in the complex-mass scheme. All real bremsstrahlung 
amplitudes have been checked by comparing different calculations for several phase-space 
points. 

2.3.2 Treatment of IR singularities 

Owing to IR singularities all the NLO pieces in (12. 2p are singular, but their sum is 
finite. When combining virtual and real corrections, soft singularities as well as singu- 
larities connected to collinear configurations in the final state cancel for "infrared-safe" 
observables after applying a jet algorithm. The remaining singularities connected to the 
collinear initial-state splittings factorize and are removed via MS PDF redefinitions, i.e. 
they are compensated by the term dafact.ab in (12. 2p . 

To isolate the IR divergences and cancel them analytically, we employ the subtraction 
formalism. To this end, we rewrite (12. 2p as 

dcTab = / d(To,ab + / [davirt.aft - da'^^^^J + / [do-rcal.afe - dcT^cXab] 

J Je Je Ji 

+ dx [do-fact,ab - da^a"i„J 
Jo Je 

= [ dao,at + [ da«"rt.a. + / da^tu. + [ dx [ da^,^,, . (2.51) 
Je Je Jr Jo Je 

Here the added subtraction terms do"™^ da^^j ^f,, and daf^^^^ „^ are constructed in such a 
way that their sum vanishes exactly and that they individually cancel all IR singularities 
of the corresponding original terms locally in phase space. As discussed in Section [2^31 
when including FwW effects, the DPA is applied to the subtracted virtual corrections, i.e. 
the second term on the r.h.s. of (I2.5ip . otherwise there would be a mismatch in the IR 
structure between virtual corrections and corresponding subtraction terms (see Ref. [85] 
for the same reasoning concerning photon radiation). 

For the subtraction terms, we choose in-house implementations of the dipole subtrac- 
tion formalism [551458] for NLO QCD calculations. Specifically this is done in dimensional 
regularization with strictly massless light quarks (including b quarks) and, alternatively. 
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in a hybrid scheme with small quark masses and massless gluons with the respective dipole 
subtraction terms from Ref . ; the typical agreement between the results from the two 
alternative treatments is at the level of 12 digits. 

Soft and coUinear singularities in the "endpoint part" (the / operator of Refs. [55|[58]) 
of the subtraction function, dalf^^ regularized using the same regularization prescrip- 

tion (dimensional or with small quark masses) as the corresponding virtual corrections. 
No regularization is needed in the subtraction terms for the real corrections, da^^^^^f,. For 
both the qq and gg channels 12 different dipole subtraction terms need to be included 
while each gg 'channel requires only 6, since we demand b quarks with finite transverse 
momentum in the final state. 

Version 1 of the real corrections employs a fully automatized implementation of dipole 
subtraction terms, which has been developed and tested during the calculation of VY+jet 
production [321 EHl ED] • The underlying matrix elements are based on helicity ampli- 
tudes employing the Weyl-van-der-Waerden formalism of Ref. [TS], which facilitates the 
treatment of helicity correlations. The colour structure of the dipole terms is generated 
automatically from the inserted colour operators and the underlying Born amplitudes. 

In version 2 all dipole subtraction terms have been implemented by hand into the 
Monte Carlo generator. Colour and helicity correlations that enter the subtraction pro- 
cedure are generated by means of AutoDipole [HI] in NwWA and OpenLoops [HS] for 
the FwW case. 

The cancellation between real matrix elements and dipole subtraction terms has been 
verified numerically in all soft and collinear regions. The individual dipole subtraction 
terms, the subtracted real matrix elements, and the integrated subtraction terms (P and 
K terms of Refs. [SJIEH]) have been compared point- wise between the two independent 
calculations. The agreement was generally at the level of 13—14 digits. 

2.3.3 Phase-space integration and numerical performance 

To achieve sufficient numerical stability, we perform the 17- dimensional phase-space 
integration using multi-channel Monte Carlo techniques [59] with adaptive weight opti- 
mization For each of the bremsstrahlung Feynman diagrams a corresponding channel 
is taken into account in the Monte Carlo integration. In addition, the integration of the 
dipole-subtracted 2 — )■ 7 contributions is optimized by means of extra channels corre- 
sponding to the dipole kinematics. For each of the 12 different dipoles and each tree-level 
diagram a new channel is introduced resulting in 12 x 31 = 372 and 12 x 14 = 168 extra 
channels in the NwWA in the gg and qq channels, respectively. When including FwW 
effects the number of additional channels grows to 12 x 79 = 948 and 12 x 38 = 456, 
respectively. In the gg 'channel dipoles based on both gg and qq Born subprocess are 
needed. So 3 x 31 -|- 3 x 14 = 135 and 3 x 79 -|- 3 x 38 = 351 channels are added here in 
the NwWA and with FwW effects, respectively. These additional channels lead to some 
improvement in the convergence of the Monte Carlo integration. 

In version 1, the phase-space integration, implemented in C+ + , is based on 
Racoon WW [00], but the phase-space mappings are built up in a more generic way 
very similar to the approach of LusiFER |61j . 
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The Monte Carlo generator of version ^ is a further development of the one used in 
COFFER77 [13] and for the calculation of the NLO corrections to pp Hjj + X 
and pp ^ttbb + X [7T] . 

For the complete NLO cross section we found agreement between the two versions 
of our code at the few-per-mille level. We also performed a detailed comparison of all 
differential distributions presented in this paper and found good agreement. 

For a typical run we generate 10® phase-space points, and the fraction of points within 
cuts is roughly 50%. This yields an accuracy for the NLO cross section of about 0.2-0.3%. 
For the LHC, the contributions of the gu^ g'd' uu, and dd channels were calculated for 
every 2nd, 5th, 4th, and 5th event, respectively. The virtual corrections for the gg, uu, 
and dd channels were evaluated for every 100th, 200th, and 500th event, respectively. 
The resulting runtime on a 3 GHz Intel Xeon processor using the if ort compiler is about 
170 h. The bulk of the runtime is taken by the real corrections to the gg channel. For the 
virtual corrections the CPU time is dominated by the gg channel and amounts to 180 ms 
per event. For the Tevatron, the weights of the different contributions are different, the 
runtime is somewhat smaller for a comparable accuracy. For the numerical results shown 
in the following almost 10^ events were generated. 

3 Numerical results 

In the following we present predictions for various observables that are relevant for top- 
pair production, either as signal or as background to Higgs production or new physics. 
Neither the following setup nor the considered observables are deliberately chosen to 
enhance the calculated off-shell and finite-top-width effects. Our aim rather is to quantify 
these effects for the standard observables. 

3.1 Input parameters 

We study the process hih2 — )■ z/ee+/i~z/^bb + X at Tevatron (pp collisions) with a/s = 
1.96 TeV and at the LHC (pp collisions) for the collider energies 7 TeV, 8 TeV, and 14 TeV. 
The input parameters and the default setup are basically the same as in Ref. |2I], where 
first results of our calculation have been presented. For completeness we specify the input 
and discuss modifications. 

In NLO (LO) QCD we employ the MSTW2008NLO (LO) parton distributions [95] and 
use the running of the strong coupling constant as with two-loop (one-loop) accuracy as 
provided by the LHApdf library. In the renormalization of the strong coupling constant 
the top-quark loop in the gluon self-energy is subtracted at zero momentum, and the 
number of active flavours is N-p = 5. In this scheme, the running of ag is generated solely 
by the contributions of the light-quark and gluon loops. Contributions induced by the 
strongly suppressed bottom-quark parton density are neglected. For the gauge-boson and 
top-quark masses we use 

mt = 172.0 GeV, Mw = 80.399 GeV, Mz = 91.1876 GeV. (3.1) 

The masses of all other quarks, including b quarks, are neglected. In view of the negligibly 
small Higgs-mass dependence we adopt the Mh 00 limit, i.e. we neglect closed fermion 
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loops involving top quarks coupled to Higgs bosons. The electroweak coupling is derived 



from the Fermi constant in the G^-scheme, 



= 1.16637 X 10"^ GeV-^ a = ^G^M^ f 1 - ^ ) • (3-2 



TT 



'Z 



While in narrow- W-width approximation (NwWA) we employ the usual real- valued weak 
mixing angle, = 1 — M-^/M|, when including finite-W-width (FwW) effects we employ 
the complex- valued mixing angle (12.281) . To derive the electromagnetic coupling (13. 2p . we 
always use the real W- and Z-boson masses. 

The leptonic decays of W bosons are treated in two different ways: in the (spin- 
correlated) NwWA and including FwW effects. In the latter discussed in 
Section I2.1.3[ FwW contributions have to be included both in the matrix elements and in 
the top-quark width used as input parameter. The top-quark width for unstable W bosons 
in NLO QCD was given in Ref. [HBj- Neglecting the bottom-quark mass, as we do through- 
out this paper, the LO top-quark width reads 



- iQV2nm^ Jo (1 - y/yy + 7^ ^'^'^ ^^'^^ 

with 7w = Fw/Mw, y = (Mw/?Tit)^, and 

Foiy)=2il-yfil + 2y). (3.4) 
Including NLO QCD corrections, the top-quark width is given by 



Gf,m^ f dy 7w 



2^ 

Fo{y) - -^F^iy) (3.5) 



with 



F,{y) = 2(1 - y)\l + 2y) [n^ + 2 Lisd/) - 2 Li2(l - y)] 

+ 4y(l -y- 2y2) ln{y) + 2(1 - yfib + Ay) ln(l - y) 

-(l-y)(5 + 9y-62/2). (3.6) 

In NwWA, i.e. for 7-w — 0, the top width follows from (I3.3l) - fl3.6l) by the replacement 

I 2 ^T^yKy-y)- (3-7) 

(1 - y/yr + iw 

With the above formulas and our input parameter set we obtain 

Tl^ = 1.4655 GeV, Ff^o = 1.3376 GeV (3.8) 

in the NwWA and 

Fj^° = 1.4426 GeV, Ff^° = 1.3167 GeV (3.9) 

including FwW corrections. Since the leptonic W-boson decays do not receive NLO QCD 

corrections and z/ee"'"yU~z/^bb production does not involve Z — )■ // subprocesses, for the 
gauge-boson widths we use the NLO QCD valued 

Fw = 2.09974 GeV, Fz = 2.50966 GeV (3.10) 

everywhere, i.e. for LO as well as for NLO matrix elements. 



* Using the measured Z-boson width instead has no significant effect on our results. 
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3.2 Jet definition, cuts, and scale choice 



We now turn to the event selection. Final-state quarks and gluons with pseudo- 
rapidity |?7| < 5 are converted into infrared-safe jets using the anti-/cT algorithm [97] . 
The jet-resolution parameter R is set to R = 0.4 and R = 0.5 for the Tevatron and the 
LHC, respectively. After recombination, we impose cuts on the transverse momenta and 
pseudo-rapidities of the leptons and b jets, and on the missing transverse momentum. For 
Tevatron we choose 

PT,b > 20 GeV, Ir^bl < 2.5, pT,miss > 25 GeV, (3.11) 

PT,z>20GeV, |r7i|<2.5, 

and for the LHC 

PT,b > 30 GeV, Ir^bl < 2.5, pT,miss > 20 GeV, (3.12) 

PT,z>20GeV, |r/d<2.5, 

where pT,miss is obtained from the vector sum of all visible transverse momenta after jet 
recombination. 

For the factorization (/ip) and renormalization (/ir) scales we have considered different 
choices. A common approach, which was adopted in Refs. [3X1132] . is to set /iR = /ip = '"^t- 
However, as shown in Section [33| this choice leads to perturbative instabilities in the 
high-energy tails of differential distributions. The use of a dynamical scale, which adapts 
to the hard scattering energy, guarantees a much better convergence of the perturbative 
expansion. We thus consider an alternative scale choice, based on the kinematic variable 



= y ^m? + p^.^m^, + (3.13) 

which corresponds to the geometric average of the top- and antitop-quark transverse 
energies. Similar scales had already been used in early papers [UIHIl] . The top transverse 
energy E^^ coincides with rrit for vanishing transverse momenta and adapts to the higher 
scattering energy at large transverse momenta. While we use the scales /igx = rrit or 
/idyn = Et to describe Uge'^ fi'u^hh production at Tevatron, for the LHC we use half of 
these scales, i.e. //fix = mt/2 or /Xdyn = -Et/2. The different scale choice is motivated by 
the fact that tt production at the Tevatron is dominated by s-channel quark-antiquark 
annihilation, while the dominant tt production mechanism at the LHC is t-channel gluon 
fusion, which prefers smaller scales. Moreover, in Ref. [98] it has been demonstrated 
that the contributions beyond NLO in an NLL soft-gluon resummation are smaller for 
fi = mt/2 than for n = rrit. Our different scale choices are summarized in Table [H There 
we also show the logarithmic average fldyn of the dynamical scale, defined via 

In^.,. = (3.14) 

The numerical values of fidyn indicate that, for what concerns the integrated cross section, 
using the dynamical scale corresponds to an effective increase of the fixed scale by roughly 
18—27%, depending on the collider energy. 



25 



poll inpy 


/ Id 


1 1 J 


^ /I \ Tp VI 


pdyn [ V J 


Ic VdLlUll 


mt 




i .yu 


zuo. i 


LHC 


mt/2 
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14 
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Table 1: Fixed (/isx) and dynamical (yUdyn) scales used for Tevatron and LHC predictions. 
The last column shows the logarithmic average of the dynamical scale, as defined in (13.141) . 



The scale uncertainty of our LO and NLO predictions is determined by uniform vari- 
ations of the renormalization and factorization scales, 

yWF = /^R = At, (3.15) 

around the central values 



fJ'O = Atfix or fiQ = /idyn- (3.16) 

When varying the renormalization scale in PDFs and matrix elements, we keep it fixed in 
the top-quark width, which is always evaluated at the scale mt. The mismatch between 
the scales used in partial and total top-decay widths is compensated by the partial-width 
correction fl2.27p . which we include in NLO^ predictions as discussed in Section [2.1.21 
In Section [231 to investigate the scale dependence of the LO and NLO integrated cross 
section we vary fi up and down by a factor eight. For all other results we provide LO and 
NLO predictions with uncertainties corresponding to factor-two scale variations. More 
precisely, the observables are evaluated at three different scales, /i//io = 0.5,1,2; the 
central value is obtained for /i = /iq, and the error band is determined by the envelope of 
the three scales. 



3.3 Total cross section and scale dependence 

We first present Tevatron and LHC total cross sections with fixed and dynamical 
scales (see Table [1]), both in NwWA and including FwW effects. The results, listed in 
Table [21 correspond to the standard cuts defined in (13.111) and (I3.12p . Besides the LO 
and NLO u^e'^ fi^U/jbh cross sections, we also present improved results that include the 
NLO"*" matching corrections of Section [2.1.21 and can be regarded as our best predictions. 
These o"^g+^-p 1,5 results incorporate higher-order effects determined by the relation (12.251) 
between the narrow-top-limit of the inclusive Uf.e'^ fi~u^hh cross section and the on-shell 
tt cross section, as well as the "partial- with correction" (I2.27p . 

At the Tevatron, the K factor is 0.95 for the fixed scale and 1.03 for the dynamical 
scale. Similarly, for the LHC we find K factors slightly below one for the fixed scale and 
5—12% above one for the dynamical scale. If we had not reduced the central scale by a 
factor 2, the K factors at the LHC would be around 1.27—1.38. The use of the dynamical 
scale increases the K factor by about 10% and changes the NLO predictions by less than 
1%. The inclusion of the NLO corrections reduces the scale dependence from 36—45% to 
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Table 2: Total cross section for pp/pp — > Vc^^ ^'v^hh + X within cuts for Tevatron and 
LHC at different CM energies in the NwWA and including FwW effects both with fixed 
and dynamical scales. In the cross-section numbers the upper variation corresponds to 
the scale = 0.5/io, the lower to = 2/io. 

4—13% at the LHC and from 44% to 9% at the Tevatron. In most cases, the NLO cross 
section goes down irrespectively of the direction of the scale variation. This indicates the 
presence of a maximum of ctnloI/^) in the vicinity of the central scale /xq (see Figures |6HZ1) . 
The residual scale dependence is comparable for fixed and dynamical scales. Including 
the partial-width correction fl2.27p in the fixed-order NLO predictions would result in a 
shift of roughly -|-2% in the central values of the LHC cross sections (due to the reduced 
central scale), while the upper (/x = 0.5/io) and lower (/i = 2/io) NLO scale variations at 
Tevatron and the LHC would move by -|-2% and —2%, respectively. 

The FwW effects amount to —0.3% at the Tevatron and —0.4% at the LHC rather 
independent of the collider energy and the scale choice. These per-mille-level FwW effects 
confirm the strong suppression anticipated in Section I2.1.3[ If we did not use the top width 
calculated with FwW effects, the inconsistency of the branching fraction (12.101) would lead 
to a fake FwW effect of roughly —3%. The NLO"*" correction, which matches our NLO 
results to the on-shell tt cross section and includes the partial- width correction (12.271) . 
ranges from —4% to -|-1%, depending on the scale choice and the collider energy. 
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Figures |6] and [7] display the dependence of the LO and NLO cross sections under 
uniform variations of the fixed and dynamical QCD scales at Tevatron and the LHC with 
a/s = 7TeV, 8TeV, and 14TeV. The fact that the central scale /iq is shifted to the left in 
the LHC plots reflects the collider-dependent scale choice (see Tabled]). For Tevatron, the 
maximum of the NLO cross section is obtained for yU ~ /xq = mt or i^T, and the K factor is 
close to one in this region. For the LHC, the maximum of the NLO cross section is shifted 
to smaller values near /x ~ /xq = ""^t/S or £'t/2, and also the K factor reaches a minimum 
there. These features are quite independent of the LHC energy and can be attributed 
to the different dominant production channels at Tevatron and the LHC. The improved 
stability of the cross section against NLO corrections and scale variations supports the 
use of the "reduced scales" /xgx = °^tl'^ and /Xdyn = -Et/2 at the LHC. 

3.4 Asymmetries 

Asymmetries in top-antitop production have found particular interest recently, trig- 
gered by experimental results from Tevatron [99] - ll01j . The situation is different for pp 
collisions at Tevatron and pp collisions at the LHC. In this section we discuss FwW effects 
and compare asymmetries obtained with fixed and dynamical scale choices. Results based 
on dynamical scales and including FwW effects should be regarded as our "best" NLO 
QCD predictions. For a recent review, including a discussion of electroweak correction 
effects [inilinS], we refer to Ref. [103] . 



S.J^.l Asymmetries jor Tevatron 

At LO, on-shell top-antitop production is totally symmetric under the exchange t f-)- 1 
both for qq and gg production mechanisms. For the gg channel this is a consequence of 
C invariance of QCD, for qq the symmetry is accidental. Consequently, the LO angu- 
lar distributions of the top and antitop quarks are forward-backward symmetric for pp 
collisions. This accidental symmetry is violated by NLO corrections, in particular by 
interference terms in qq — )■ ttg and gg — ?■ ttg radiative processes |104j and by interfer- 
ence terms in the qq channel between virtual corrections and the leading order [105j. 
Since top-antitop production at Tevatron is dominated by the qq channel, an observable 
forward-backward asymmetry emerges. In the laboratory frame this can be defined as 

.FB _ (y{yt > 0) - (y{yt < 0) .^lyx 
- a{y,>0) + aiy,<Oy ^^'^^^ 

where yt and yi are the rapidities of the top and antitop quark, respectively. In our 
predictions we set a = cjnlo both in numerator and denominator. Owing to CP symmetry, 
the asymmetry of the antitop is oppositely equal to the asymmetry of the top quark. 
Our numerical results for the asymmetry f l3.17p are given in the second column of the 
upper part of Table [3] including FwW effects. Compared to the result for the fixed scale 
fiQ = rrit, the asymmetry is lowered by 10—15% for the dynamical scale /xq = Et- The 
scale uncertainty is large, since the numerator is effectively known only in leading order 
and has a scale dependence like a^, while the NLO cross section in the denominator is 
rather stable against scale variations. Within the integration accuracy, the results for 
on-shell and off-shell W bosons coincide. 
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Figure 6: Scale dependence of the LO and NLO cross sections at Tevatron (left) and the 
LHC with ^/s = 7TeV (right). The renormalization and factorization scales are varied 
around the fixed (/xq = /^fix) or dynamical (/xq = /^dyn) central values defined in Table [H 




Figure 7: Scale dependence of the LO and NLO cross sections at the LHC with ^/s = 8 TeV 
(left) and 14 TeV (right). The renormalization and factorization scales are varied around 
the fixed or dynamical central values /xq defined in Table [H 
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In the literature, the asymmetry A™ab often been calculated by using the LO 
cross section in the denominator pi H [29 |[TU2lllU5] . This reduces the scale dependence by 
at least a factor of 4, while using NLO predictions in the numerator and the denominator 
of (13.171) . as in the present paper, provides a more conservative estimate of theoretical 
uncertainties. In any case, this choice does not really affect the size of the finite-width 
effects that we are after. We do not include QED or electroweak corrections. 

We note that for off-shell top quarks an asymmetry already exists in LO, e.g. from non- 
resonant diagrams where W bosons are emitted from initial-state quarks. This asymmetry 
is, however, only at the level of 10~^ and thus negligible. 

An analogous asymmetry can be defined for the charged leptons. Results for the 
positron including FwW effects are given in the third column of the upper part of Table [3l 
both for fixed and dynamical scales. Our results at the fixed scale /io = rrit for both lab 
and LAB ^grse well with the corresponding numbers given in Refs. P ^ [52 |[TUB] . 

Alternatively, asymmetries can be defined in the CM frame of the tt system as 

.FB _ crjyt > yi) -(y{yt< yi) 

^•"^ " a{y.>yi) + a{y,<yiy ^^^ ^^ 

At the CDF experiment, these asymmetries have been measured separately in two regions 
of the tt invariant mass Mti, viz. for Mtt > 450 GeV and M^t < 450 GeV, with the 
result [TOT] 

= 0.162(47), 

A™M°^(^tt < 450 GeV) = 0.078(54), 

A™jS°^(Mtt > 450 GeV) = 0.296(67). (3.19) 
The DO experiment has measured |100] 

= 0.196(65). (3.20) 

From our calculation, using /io = -Et and including FwW effects, we find the results given 
in the lower part of Table [31 While these results agree well with theoretical calculations 
of other groups [2 ^1106] , the measured asymmetry is significantly higher, in particular for 
Mtt > 450 GeV. 



3.4-2 Asymmetries for the LHC 

At the LHC, the forward-backward asymmetry of the quark- and antiquark-induced 
partonic subprocesses does not show up in observables, since all distributions are forward- 
backward symmetric for the proton-proton initial state. However, as a result of the 
dominance of valence quarks over sea (anti) quarks at large momentum fractions, the 
partonic asymmetry manifests itself as a hadronic central-edge asymmetry. In practice, 
in the laboratory frame, antitop quarks tend to be produced more centrally than top 
quarks. The ATLAS collaboration [107] measures an asymmetry based on the LAB-frame 
rapidity of the top quarks. 
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,, aFB aFB 

^t,LAB ^e+,LAB 

^+0.034 n nQRI /'c;^+0.026 



hO.019 
-0.009 



Et 0.0454(5)1°:°26 0.0321(5)+°-°^9 



/xq /L^^m A^gM(^tt < 450GeV) A%(^tt > 450 GeV) 

m, 0.0749(5)1°:°^} 0.0491(6)1°;°?^ 0.1281(11)1°:^! 
0.0683(5)1°:°?^ 0.0486(6)1°:°?^ 0.1078(9)1°:°^^ 



Table 3: Forward-backward asymmetries for the top quark and the charged lepton at the 
Tevatron including FwW effects, using a dynamical /iq = or a fixed = rrit scale. 
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14 
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mt/2 
Et/2 
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Table 4: Central-edge asymmetries for the top quark and the charged lepton for the LHC 
at different CM energies including FwW effects, using a dynamical /iq = -Et/2 or a fixed 
/io = mt/2 scale. 



while CMS [108J prefers to investigate an asymmetry based on the pseudo-rapidity, 

^CE,,, ^ ^(l^tl > \Vt\) - (^{\Vt\ < \Vt\) ^2 22) 



Again these asymmetries can be defined for charged leptons as well, 

^CE,y ^ cr{\yc+\ > |y^-|) -(t(|i/c+| < \y^-\) 



Since the charged leptons can be considered massless, Af^'^ = A'^^'^. Results for the 
top-quark and charged- lepton central-edge asymmetries for the dynamical /iq = Et/2 
and fixed Hq = 'mt/2 scale and with FwW effects are given in Table |H Also in this case 
FwW effects turn out to be strongly suppressed (at the sub-per-mille level). 



3.5 Differential distributions 



In the following we present various distributions obtained by applying the cuts spec- 
ified in Section 13.21 For each of the observables illustrated in Figures [81430] we present 
three plots. The left plots display absolute LO (blue, dashed) and NLO (red, solid) pre- 
dictions together with corresponding uncertainty bands resulting from scale variations 
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within 0.5 < fi/fiQ < 2. In the upper-right plots we show the same LO and NLO pre- 
dictions normahzed to LO results at the central scale, i.e. Ki^o{^) = d(jLo(/^)/do"Lo(Ato) 
and Kf^i^o{fi) = dcrNLo(/^)/do"Lo(Aio)- Here the blue band illustrates the relative scale 
uncertainty of the LO cross section, and the central curve of the red band corresponds 
to the usual NLO correction factor, K = A'nlo(/^o)- These first two plots always include 
FwW effects, whose impact is illustrated in the lower-right plots. There we display the 
FwW correction factor 

^F-w = ^^NwWA - 1' (3-24) 

obtained by comparing the NwWA and the FwW variants of our calculation. 

All results in Figures [T0] - [30] are obtained with the dynamical scales /io = -Et/2 at the 
LHC and /xq = -E't at Tevatron. To motivate this choice, we first show, in Figures [8] and 
[HI that using a fixed scale = mt/2 at the LHC leads to serious perturbative instabilities 
in the high-energy tails of distributions. 

3. 5. 1 Differential distributions for the LHC at 8 TeV 

We first provide results for the LHC at 8 TeV. As examples for typical px distributions 
we present the distributions in the transverse momenta of the positron and of the top 
quark for the fixed scale /x = mt/2 in Figures [H] and [S The tails of these distributions 
are relevant for new-physics searches based on boosted top quarks, while the lepton-px 
distribution plays an important role for the tt cross-section acceptance. The typical lepton 
Pt is below 100 GeV and the average of the top quark is around 100 GeV. The FwW 
effects are within integration errors. In the plotted range, the cross section falls by almost 
four orders of magnitude, while the K factor drops by 75%. In the high-px tails, the 
corrections become so large that NLO predictions move outside the LO band and their 
scale variation exceeds —100%. This pathologic behaviour at large px can (at least in 
part) be attributed to large logarithms associated with the running of Og, which can be 
effectively resummed by adapting the QCD renormalization scale to the characteristic 
scattering energy. This motivates us to adopt the dynamical scale /io = -Et/2 as our 
standard LHC scale choice, which is applied to all LHC observables presented in the 
following. 

The PT,e+ and pT,t distributions, displayed in Figures [ID] and [HI show that the dy- 
namical scale choice leads to a drastic improvement of the perturbative stability. Now 
the K factor changes only slightly with px (within 20%), and the NLO band lies within 
the LO band. 

Next we consider the distributions in the transverse momenta of the harder b jet 
(the one with larger px) and the softer b jet in Figures [T2] and [T31 The corrections to 
the transverse-momentum distribution of the harder b jet are small near the maximum 
of the distribution and increase slowly for higher px by about 40%, while those for the 
distribution of the softer b jet decrease by 20%. Using the fixed scale /iq = mt/2 (not 
shown), the K factor to the px,b,max distribution decreases by 40% and the one for the 
PT,b,min distribution even by 90%. FwW effects are at the level of 1% or below. 

The distribution in the transverse momentum of the bb pair is presented in Figure [TH 
This variable plays an important role in boosted-Higgs searches with a large tt background. 
In particular, the strategy proposed in Refs. |109[I110] to extract a pp — )■ H(— >■ bb)W 
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Figure 8: Transverse-momentum distribution of the positron with standard cuts for the 
LHC at y/s = 8 TeV for fixed scale /iq = mt/2. 
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Figure 9: Transverse-momentum distribution of the top quark with standard cuts for the 
LHC at a/s = 8 TeV for fixed scale /xq = mt/2. 
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Figure 10: Transverse-momentum distribution of the positron with standard cuts for the 
LHC at y/s = 8 TeV for dynamical scale /io = Et/2. 
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Figure 11: Transverse- momentum distribution of the top quark with standard cuts for 
the LHC at ^/s = 8 TeV for dynamical scale /iq = Et/2. 
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Figure 12: Transverse-momentum distribution of the harder b jet with standard cuts for 
the LHC at ^/s = 8TeV for dynamical scale /iq = Et/2. 
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Figure 13: Transverse-momentum distribution of the softer b jet with standard cuts for 
the LHC at ^/s = 8 TeV for dynamical scale fiQ = 
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Figure 14: Transverse-momentum distribution of the bb pair with standard cuts for the 
LHC at y/s = 8 TeV for dynamical scale /io = Et/2. 




Figure 15: Missing-transverse- momentum distribution with standard cuts for the LHC at 
y/s = 8 TeV for dynamical scale /xq = -E't/2. 
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signal at the LHC is based on the selection of boosted H — > bb candidates with PT.bB > 
200 GeV, which permits to reduce tt contamination (and other backgrounds) in a very 
efficient way. As can be seen from Figure [T31 the suppression of tt production is indeed 
particularly strong at PT,bb ~ 150 GeV. This is due to the fact that, in NtWA, b quarks 
need to be boosted via the pt of their parent (anti)top quarks in order to acquire pT,h > 
[m^ — M^)/(2mt) ~ 65 GeV, and a bb system with high is kinematically strongly 
disfavoured at LO, since top and antitop quarks have opposite transverse momenta. The 
NLO corrections undergo less stringent kinematic restrictions, resulting in a significant 
enhancement of WWbb events at large p^hb- This is clearly reflected in the differences 
between the LO and NLO curves in the left plot of Fig. [TH At NLO the tt system can 
acquire large transverse momentum by recoiling against extra jet radiation. As indicated 
by the upper-right plot, the NLO correction represents 50—80% of the cross section at 
high Pt, corresponding to a huge i^-factor of 2—4. FwW effects (lower-right plot) stay 
at the level of 2%. On the other hand, FtW effects (not shown here explicitly) become 
as large as 10-30% for PT,bb > 200 GeV |62] . This is most likely due to non-resonant 
topologies with direct bb production from a high-px gluon that recoils against a W"^W~ 
system and splits into a bb pair. 

Figure US] displays the distribution in the missing transverse momentum, pT,miss = 
Ipt,!/c + PT,pp I • The tail of this distribution, which is relevant for new-physics searches 
based on missing transverse energy plus jets and leptons, features a qualitatively similar 
behaviour as in the case of PT,bb5 owing to analogous kinematic constraints. However, 
in the case of PT.miss the corrections are less pronounced: the K factor only rises to 2.5. 
FwW contributions are not significant. 

The distribution in the invariant mass of the top quark, = M^^^^+y^, in the vicinity 
of its resonance is shown in Figure [161 Below the resonance the NLO corrections develop 
a radiative tail of positive corrections due to final-state gluon radiation that is not recom- 
bined with the top-quark decay products. The NLO corrections induce a small shift in the 
peak of the distribution, whose magnitude depends on the jet recombination. Note that 
the major part of this shift would also be accounted for in a pure parton-shower approach 
to describe higher-order QCD effects. Thus, a thorough assessment of the impact of the 
full NLO QCD calculation with off-shell top quarks on observables that are used in the 
top mass measurement requires a careful comparison with predictions based on parton 
showers. Note also that FwW corrections increase in the side bands of the resonance to 
the level of a few per cent. 

Figure [TT] displays the distribution in the invariant mass of the positron and a b jet, i.e. 
the visible products of a top-quark decay. We here use the Monte Carlo truth to select the 
negatively charged b-quark jetl^ In narrow-top-width and LO approximation this kine- 
matic quantity is characterized by a sharp upper bound, M^f^ < — M^ ~ (152 GeV)^, 
which renders it very sensitive to the top-quark mass. The value of rrit can be extracted 
with high precision using, for instance, the invariant-mass distribution of a positron and 
a J/ip from a B-meson decay |lllllll2j . an observable that is closely related to Me+b. In 
the region below the kinematic bound, the NLO corrections to Mg+b vary between —10% 
and 20%, and the impact of the NLO shape distortion on a precision rrit measurement is 
certainly significant. For Me+b < 150 GeV, the NwWA is very good. Above the kinematic 

^In Ref. [28] the e^b pair is formed by selecting the b jet that yields the smallest invariant mass. 
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Figure 16: Invariant-mass distribution of the top quark, Mt = M,^^e+b5 with standard cuts 
for the LHC at ^/s = 8TeV for dynamical scale /io = -E't/2. 




Figure 17: Invariant-mass distribution of positron-b-jet system with standard cuts for 
the LHC at ^/s = STeV for dynamical scale /io = -Et/2. 
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Figure 18: Distribution in the total transverse energy Ht with standard cuts for the LHC 
at a/s = 8 TeV for dynamical scale /io = Et/2. 



bound, NLO corrections become clearly visible, giving rise to a tail that extends above 
^e+b ~ "^t ~ ^^"^ ^^^'^ FwW corrections become sizeable. In this kinematic region 
also the finite top width causes effects at the level of 50% [62]. While the contribution to 
the total cross section from the region above 150 GeV is fairly small, the impact of these 
contributions on the top-mass measurement might be non-negligible, given the high rrit 
sensitivity of the M^+^ ~ — region. Again a careful comparison between NLO 
off-shell calculation and parton-shower approach would be required to quantify off-shell 
effects on the mt measurement. 

In Figure HE] we display the distribution in the total transverse energy, 

Ht = PT,e+ + PT,m- + PT.miss + X^Ptj, (3.25) 

j 

where the scalar sum over the jet transverse momenta, Ptj, includes the two b-quark jets 
as well as possible extra jets contributing to NLO real emission. The Ht distribution is 
relevant for new-physics searches. The large correction above Ht = 400 GeV results from 
the inclusion of the gluon jet, which shifts the distribution to higher Ht- If this is not 
included, as e.g. in Ref. [32], the K factor is flat and varies between 0.8 and 1.2. 

Finally, we show the distributions in the invariant masses of the bb and tt pairs, in 
Figures [19] and [201 and the distributions in the azimuthal angle between the two charged 
leptons in the transverse plane and the cosine of the angle between them, in Figures [51] 
and [22] In general these distributions receive moderate shape corrections, and FwW 
effects are tiny. Sizable NLO and FwW effects show up only in the Mtt distribution below 
the tt threshold. 
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Figure 19: Invariant-mass distribution of the bb pair with standard cuts for the LHC at 
^/s = 8TeV for dynamical scale = Et/2. 
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Figure 20: Distribution in the invariant mass of the tt pair with standard cuts for the 
LHC at ^/s = 8TeV for dynamical scale /iq = -E't/2. 
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Figure 21: Distribution in the azimuthal angle between the positron and the muon in the 
transverse plane, 4>e+fj.-, with standard cuts for the LHC at ^/s = 8TeV for dynamical 
scale /io = Et/2. 
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Figure 22: Distribution in the cosine of the angle between the positron and the muon, 
cos9e+fj,-, with standard cuts for the LHC at y/s = 8TeV for dynamical scale /iq = Et/2. 
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3.5.2 Differential distributions for the LHC at 14TeV 



At the LHC with 14TeV, using the dynamical scale /io = -Et/2, we observe that 
the shapes of the distributions and the NLO corrections look qualitatively similar as for 
8 TeV. However, the K factors for transverse-momentum distributions grow faster with 
Pt than at 8 TeV. In Figures |23] and [2l] we show results for the distributions in the 
transverse momenta of the positron and of the harder b jet. Above px = 100 GeV, the 
K factor grows by about 30% for the PT,e+ distribution and by about 60% for the PT,b,max 
distribution. Using the fixed scale = mt/2 instead (not shown), the K factor decreases 
by about 40% for the PT,e+ distribution and by about 20% for the PT,b,max distribution. 
In general, FwW corrections remain similarly suppressed as at 8 TeV. 

3. 5. 3 Differential distributions for Tevatron 

Next we show a few distributions for Tevatron, using the dynamical scale /ig = -Et- 
The transverse-momentum distribution of the positron is presented in Figure [221 The K 
factor decreases slightly with increasing p^. FwW effects are completely negligible. 

We turn to the distributions in the transverse momenta of the harder b jet and the 
softer b jet in Figures [26] and 1271 Below the maxima of the distributions we find a positive 
K factor at the level of 40%, but above pT,b = 70 GeV the K factor is pretty flat and 
close to one. For large PT,b,min FwW effects at the level of 2% show up. 

The distribution in the transverse momentum of the bb pair is presented in Figure [2S1 
In contrast to the LHC there is no enhancement of WWbb events at large Px.bB; cind the 
K factor stays near one in the complete considered range. A similar behaviour can be 
observed for the missing-transverse- momentum distribution (not shown). 

The distribution in the invariant mass of the top quark in the vicinity of its resonance 
is shown in Figure [23 While the enhancement of the distribution by NLO corrections 
below the resonance is similar as for the LHC, the rise of the K factor above the resonance 
is weaker. In the side bands of the resonance FwW corrections are at the level of a few 
per cent as at the LHC. 

The distribution in the total transverse energy as defined in fl3.25p is presented in 
Figure [3D1 The effect of the NLO corrections is qualitatively similar as for the LHC, but 
the enhancement of the distribution for high Ht is only 40%. As for the LHC, the K 
factor would become much flatter if the gluon px was not included in H^. 

3.6 Limit of on-shell top quarks 

To quantify non- resonant and off-shell-top contributions to the z/ee^/i~z/^bb integrated 
cross section, we have investigated its narrow-top- width limit. Ft — 0, by means of a 
numerical extrapolation. This is shown in Figure [HH where we plot 



in the range < Ft < Ff where F^ is the physical top-quark width. The fac- 
tor (Ft/Fli''^-^'')^ compensates deviations of Ft from T^^^^ in such a way that the effective 




(3.26) 
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Figure 23: Transverse-momentum distribution of the positron with standard cuts for the 
LHC at y/s = 14TeV for dynamical scale /xq = -Et/2. 
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Figure 24: Transverse-momentum distribution of the harder b jet with standard cuts for 
the LHC at y/s = 14TeV for dynamical scale /iq = -E't/2. 
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Figure 25: Transverse-momentum distribution of the positron with standard cuts for 
Tevatron for dynamical scale fio = E^. 

top-decay branching fractions remain constant. The consistency of the Ft — > extrap- 
olation at LO is demonstrated by the green squares in Figure EH which were obtained 
with an explicit Born calculation in NtWA. At NLO, when Ft — )■ the virtual and real 
parts of NLO corrections are individually enhanced by soft-gluon logarithmic singularities 
oc OglnFt, which are not reshuffled via Catani-Seymour dipoles. The precise numerical 
cancellation of such singularities in the sum of virtual and real corrections manifests it- 
self in the quality of the linear convergence of the Ft — )■ extrapolation. This provides 
a non-trivial confirmation of the consistency and numerical stability of the calculation. 
Finite-top-width effects are extracted by comparing results at Ft = T^^^^ and Ft — )■ 0. At 
the Tevatron, we find that FtW effects shift the LO(NLO) cross section by -0.8%(-0.9%). 
This is fairly close to the numerical value of Ft/mt, which represents the expected order 
of magnitude. At the LHC, with ^/s = 7 TeV, FtW effects turn out to be even smaller: 
+0.4% at LO and +0.2% (comparable to the Monte Carlo statistical error) at NLO. 
Most likely this difference to the result of ~ 1% at Tevatron is due to the influence of 
singly-top-resonant diagrams corresponding to associated tW production, whose positive 
contribution even overcompensates the negative effects from off-shell tt production. 

The results of this extrapolation have been conflrmed in a dedicated comparison [62] 
between our calculation and the one of Ref. [28] for on-shell top quarks. A similar Ft — )■ 
extrapolation has been presented in Ref. |32j . 

3.7 Comparison with HELAC-NLO 

Our predictions for the hadronic i/eS^/i^t'^bb cross section have been successfully re- 
produced by the HELAC-NLO collaboration and vice versa. On the one hand, our 
results for the Tevatron in NwWA [31] have been checked by the HELAC-NLO group 
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Figure 26: Transverse-momentum distribution of the harder b jet with standard cuts for 
Tevatron for dynamical scale /xq = -^'T- 
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Figure 27: Transverse-momentum distribution of the softer b jet with standard cuts for 
Tevatron for dynamical scale /iq = -Et- 
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Figure 28: Transverse-momentum distribution of the bb pair with standard cuts for Teva- 
tron for dynamical scale /io = Et- 
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Figure 29: Invariant-mass distribution of the top quark, Mt = Mj^^g+b, with standard cuts 
for Tevatron for dynamical scale /xq = -Et- 
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Figure 30: Distribution in the total transverse energy with standard cuts for Tevatron for 
dynamical scale fio = E^. 



pp -;• Ure+fi-j^f^hh+X @ = 1.96 TeV 




pp -> Uce+ fi^ u^,hh+X @ = 7 TeV 



45 
44.5 

44 
43.5 

43 
42.5 

42 
41.5 



0.2 0.4 0.6 0.8 1 

Tt [GeV] 



LO 
NLO 




1.2 1.4 




650 



0.2 0.4 0.6 0.8 1 1.2 1.4 

Ft [GeV] 



Figure 31: Numerical narrow-top- width extrapolation of the LO and NLO z/ce+yU^P^bb 
cross section at Tevatron and the LHC at y/s = 7 TeV. This study does not include 
finite- W- width effects and is based on the same fixed-scale choice, fio = rrit, for Tevatron 
and LHC. 
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collider 


energy 


aLo[fb] 


aLo[fb] 


crNLo[fb] 


0-NLo[fb] 






Ref. [32J 


this work 


Ref. [32] 


this work 


Tevatron 


1.96 TeV 


34.922(14) 


34.923(4) 


35.705(47) 


35.777(24) 


LHC 


7TeV 


550.54(18) 


550.30(6) 


808.46(98) 


807.56(52) 


LHC 


10 TeV 


1394.72(75) 


1394.51(14) 


1993.3(2.5) 


1992.1(1.4) 



Table 5: Cross sections for off-shell top-pair production for various colliders in the setup 
of Ref. [32j with anti-fcx algorithm. 

finding agreement within integration errors |113j . On the other hand, we compared the 
FwW variant of our calculation against the HELAC-NLO predictions of Ref. [32]. The 
corresponding results obtained with the anti-/cT jet algorithm and the input parameters, 
cuts, and PDFs of Ref. [32] are listed in Table [51 At LO and NLO both calculations co- 
incide within (1— 2)cr, which is for NLO at the level of 0.2%. We note, however, that the 
numerical values of the input parameters rj"*"* and rf""° used in Ref. [32], and adopted 
for the technical comparison in Table [5l do not include FwW corrections. As discussed 
in Sections 12. 1.31 and 13. 3[ this induces a fake FwW shift of roughly —3% in the integrated 
cross section. 

4 Conclusions 

The precise understanding of top-antitop production is important for precision tests of 
the Standard Model and in searches for the Higgs boson and new physics. In this paper, 
we have presented NLO QCD predictions for the process pp/pp — Uce^ fi^u^jbh + X, 
which corresponds to tt production and decay in the di-lepton channel. The calculation 
includes all finite-top-width effects, i.e. off-shell top contributions, diagrams involving 
less than two top resonances, and all factorizable and non-factorizable NLO corrections. 
The leptonic W-boson decays have been described alternatively within a narrow- W-width 
approximation or including finite- W-width effects. 

The employed techniques, based on Feynman diagrams and tensor integrals, provide 
high numerical stability and CPU efficiency. Each Feynman diagram is algebraically 
reduced to a canonical form and automatically processed to Fortran output, while 
tensor loop integrals are reduced to scalar integrals using numerical algorithms that avoid 
instabilities due to inverse Gram determinants and other spurious singularities. The real 
corrections are integrated using the dipole subtraction method. 

We have presented a detailed study of integrated and differential cross sections and of 
asymmetries at Tevatron and the LHC at -y/i = 7 TeV, 8 TeV, and 14 TeV. We pointed 
out that a fixed QCD scale leads to perturbative instabilities in the tails of transverse- 
momentum distributions, while using a dynamical scale related to the transverse energy of 
the top quarks yields stable NLO predictions. For the LHC, where top-antitop production 
mainly proceeds via gluon-gluon annihilation, we decided to reduce the QCD scale by a 
factor two as compared to the Tevatron, where the quark-antiquark channel dominates. 
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This can be motivated by resummation of higher-order corrections, and leads to smaller 
K factors and a somewhat reduced residual scale dependence. 

We also discussed the consequences of truncating the perturbative expansion at NLO 
in presence of unstable intermediate particles. Although the NLO corrections to the pro- 
duction and decay of top quarks are separately well known, a fixed-order NLO description 
of the entire process, i.e. production and decay, does not include the product of these two 
corrections. Such contributions are of NNLO kind but can largely exceed one per cent. 
To take them into account, we introduced an approach based on the idea of matching the 
inclusive z/ee"^/i~z/^bb cross section, in the Ft — ?■ limit, to the on-shell tt cross section. 
Numerically, the corresponding correction amounts to a few per cent. 

Effects associated to the finite top-quark width are formally suppressed by Tt/rrit 
and turn out to be below 1% for the total cross sections at Tevatron and the LHC at 
a/s = 7TeV. In contrast, finite-top-width effects can reach tens of per cent in phase space 
regions where top-antitop production is suppressed as a signal, but remains important 
as a background to Higgs and new-physics searches, as discussed in Ref. |62] where our 
results are compared to the ones of the (spin-correlated) narrow-width approximation of 
Ref. [28]. 

In spite of the larger numerical value of Fw/Mw as compared to Ft/mt, finite- W- width 
effects turn out to be less important than finite-top-width effects. For observables that are 
dominated by on-shell top quarks they remain below 0.5%, and for all other observables 
considered in the present paper they never exceed a few per cent. The only exception 
is the Me+b distribution in the vicinity of M'^+^ = — M^, an observable which is, 
however, important for the precision measurement of the top-quark mass. The small size 
of finite- W-width contributions (to inclusive observables) is due to a double-suppression 
mechanism, which involves subtle cancellations between finite- W-width corrections to 
u^e'^ fi~i?^hh matrix elements and to the Ft input parameter. Besides these observations, 
one should keep in mind that observables where top-quark and/or W-boson resonances 
are kinematically excluded can receive much larger finite-width contributions. 

The step from on-shell top quarks to the inclusion of the top-quark decays with its 
full off-shell effects and the differential information on the decay products, including NLO 
QCD corrections, is an important step towards pushing the simulation of tt production 
at hadron colliders to a higher level. A natural next development, apart from extending 
the calculation to hadronic W decays as well, consists in a proper matching of the new 
NLO QCD corrections to off-shell tt production with QCD parton showers. Further 
future improvements of the calculation should also aim at including NLO electroweak 
corrections, a task that actually leads to a horrible NLO calculation for a 2 — )■ 6 particle 
process. However, our setup of using a double-pole approximation for the W resonances 
effectively factorizes the full process again into a 2 — )■ 4 (WWbb) production and the 
1 — > 2 W-boson decays, rendering this calculation possible, though still very demanding. 
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particle 


p°[GeV] 


[GeV] 




p3[GeV] 


a 


250 








250 


h 


250 








-250 




48.56688853967434 


-5.434957739109776 


-19.57544900121795 


44.11355452919156 


e+ 


44.98131949428399 


27.80765763034380 


-19.96763010585797 


-29.17785167753225 




81.94400052131301 


35.50973810316283 


57.00882135433334 


-46.94541520821219 




75.44248660279438 


-13.21633653678938 


-1.177064626981959 


-74.26649144952955 


b 


156.4430715983779 


-88.78355499018868 


-95.32442077834540 


86.63238319307789 


b 


92.62223324355639 


44.11745353258120 


79.03574315806995 


19.64382061300453 



Table 6: Phase-space point for the partonic process ah — )■ z/oe^/i~z/^bb. All scattering 
particles are massless, intermediate top quarks are (slightly) off shell, and intermediate 
W bosons are exactly on shell, i.e. {pi,^ +Pe+)^ = {p^l- +Pp,J^ = ^w- 

A Benchmark numbers for the virtual corrections 

In order to facilitate a comparison to our calculation, we provide explicit numbers 
for the squared tree amplitude and the corresponding virtual correction for the partonic 
processes uu/gg — )■ VeQ'^ l^'v^jbh evaluated at the phase-space point specified in Table El 
The amplitudes are computed using the input parameters of Section 13.11 and the complex 
top-quark mass (12. 3p . For the top-quark width we use the NLO value in (13. 9p . which 
includes FwW corrections, both in the tree and one-loop matrix elements. Since we 
compute virtual matrix elements handling W-boson resonances in double-pole approxi- 
mation, we consider a phase-space point with on-shell intermediate W bosons. Moreover, 
as discussed in Section [2.1.31 we treat the W- and Z-boson masses as real parameters, 
we set Fz = 0, and we insert a non-zero W-boson width only in the W-boson propaga- 
tor, which contributes a factor l/(FwMw) for each W resonance at = M^. For the 
renormalization scale and the scale of dimensional regularization we choose n = fi^ = rrit. 
The corresponding values of the strong coupling constant in the renormalization scheme 
described in Section [3] are 

tts(mt)|LO = 0.1258086856923967, as(mt)|NLO = 0.1095186120831399. (A.l) 

In the following we use the NLO strong coupling everywhere, i.e. also in tree matrix 
elements. 

For the spin- and colour-averaged squared LO amplitudes in DPA we obtain 

l-^uSoPAr = 1.568863069202805- 10"^ GeV"^ 

I-^ssVaP = 4.554053154627972 ■ 10"^ GeV-^ (A.2) 

while including doubly- and singly- W-resonant diagrams in the complex-mass scheme 
(I2.28P with input parameters (I3.10p yields 

l-A^uSfuuP = 1.567938932324559 • 10^^ GeV~^ 
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l-^gg^fuuT = 4.551376909529072 • 10"^ GeV~®. (A.3) 

We express virtual NLO contributions in the 2 — )■ 6 phase space as Laurent series in 
6=(4-D)/2, 

I-MdpaP = (l + cri2^l,%eA \M]fp^\\ (A.4) 

\ fc=0 / 

where we factor out the LO term in DPA and the normahzation factor 

_ (47r)T(l + e)r^(l-6) _ (AnY 

TiT^) -f(r^+ 

= (47r)T(l + e)-^e2 + 0(e3). (A.5) 





We spht the result into the two parts, 

which correspond to the contributions of renormalized loop diagrams (loops) and the / 
operator of the dipole subtraction function as defined in Ref. |58]. Here we work with 
purely massless light fermions, i.e. in pure dimensional regularization. In ^l^^pg we include 
counterterm contributions from field, coupling, and mass renormalization. The numbers 
in Table [7] have been obtained in the 't Hooft-Feynman gauge using the 't Hooft-Veltman 
variant of dimensional regularization (four-dimensional external partons). 

The agreement between our two independent versions of the virtual corrections, as can 
be seen from Tabled is typically about 11 digits at non-exceptional phase-space points. 
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